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Abstract 

Measurements are reported of the yield and spectrum of fluorescence, excited by 
a 28.5 GeV electron beam, in air at a range of pressures of interest to ultra-high 
energy cosmic ray detectors. The wavelength range was 300 - 420 nm. System cali- 
bration has been performed using Rayleigh scattering of a nitrogen laser beam. In 
atmospheric pressure dry air at 304 K the yield is 20.8 it 1.6 photons per MeV. 
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1 Introduction 



We report on a measurement of the fluorescent yield of air at wavelengths 
and pressures of interest to large area cosmic ray shower detectors. Telescope 
arrays, imaging fluorescence in large volumes of the atmosphere, continue to 
probe the spectrum of ultra-high energy cosmic rays (UHECR). The present 
work is a step in establishing the foundation for the energy scale calibration 
of UHECR observations. 

The steeply falling spectrum is very difficult to measure above ~ 10^^ eV 
(1.6 Joules per particle) because of the very low rates and our technical in- 
ability to produce and study particles of such energies in the laboratory. It 
is of great interest [1], however, because conventional mechanisms do not ex- 
plain the occurrence of such high energy particles, and this holds out the 
possibility of uncovering new physics. One possibility is acceleration by very 
energetic sources [2J. Alternatively the decay of primordial super-heavy parti- 
cles [3] might be responsible. Distinguishing between the mechanisms requires 
continued statistical and systematic improvement in the measurement of the 
UHECR spectrum. 

Recent experimental interest was focused on evidence for the Greisen-Zatsepin- 
Kuzmin (GZK) |1] cutoff in the spectrum, at about 10^° eV, expected from 
interactions with the cosmic microwave background. The large ground level 
scintillation array AG AS A [S] reported an unbroken spectrum above the limit. 
The atmospheric fluorescence array HiRes [6] , however, reported a steepening 
of the spectrum consistent with the cutoff. This was recently confirmed by 
preliminary results from the nearly completed Auger detector [7j, whose en- 
ergy scale is also set by fiuorescence measurements. The latter two arrays find 
that the spectral slope becomes shallower as energy increases above ~ 5 x 10^^ 
eV, before a final steepening. But they do not yet agree on the fiux, or the 
exact energy thresholds and slope exponents of the segments of the spectrum. 
Uncertainty in the knowledge of air fiuorescence from cosmic ray showers may 
contribute to the differences. 

Building on previous experiments, increasingly precise measurements are an- 
ticipated for the spectral shape and shower profiles above ~ 10^^ eV. This 
promises to afford a rich field of study. Inferences about the relative contri- 
butions of galactic and extra-galactic fiuxes, the cosmic distance distribution 
of sources, and UHECR composition, will become accessible [8]|9][T0] . New ex- 
periments, at various stages of development, are aimed at enhanced statistics 
with concentrated study of systematic effects [TT] . All of these include at least 
a fluorescence measurement system for atmospheric showers, and depend on 
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an improved knowledge of the atmospheric fluorescence processes. 

Fluorescence detectors can measure the cosmic ray energy by measuring the 
total light yield of the shower that it generates in the atmosphere. From the 
initial interaction high in the atmosphere, the secondaries interact again, still 
at high altitude, and as the air becomes denser, a cascade quickly builds up, 
driven by the decay of the copious neutral pions produced. The decay of the 
pions leads to the electromagnetic shower that, in fact, is the dominant vehicle 
for the absorption in the atmosphere of the initial energy of the cosmic ray. 

It is the electrons and positrons of the shower that transfer energy to the air 
molecules and cause them to emit fluorescence light. The energy transfer is well 
described by the Bethe-Bloch equation [12], and a large fraction comes from 
particles in the energy range of tens of MeV. Within the gas molecules the en- 
ergy flow is very complicated because of the number of energy levels available, 
and the pressure-dependent competition between fluorescent emission and en- 
ergy transfer in molecular collisions. Work towards the theoretical modeling 
and calculation of the fluorescent light yield is progressing p^l4|15p6] . 




The UHECR atmospheric fluorescence detectors use spherical mirrors to im- 
age the profile of the light on to arrays of photomultiplier tubes piTT] . The 
light's spectral range is restricted by filters to an ultraviolet window with low 
sky light background between 300 and 400 nm. The fluorescence light in this 
wavelength range is dominated by nitrogen emission lines, with major bands 
near 315, 337, 355, 380 and 391 nm, (95% of the light) and a few others of 
lesser intensity [T31IT7] . It is important to note that modern detectors can mea- 
sure light from very high energy showers at distances exceeding 30 km. But at 
such distances the Rayleigh scattering by air molecules becomes important, 
and, even in the limited 300 to 400 nm range, its wavelength dependence 
preferentially reduces the detection efficiency of the short wavelengths. At 337 
nm, for example, the exponential scattering length at one atmosphere is 11 
km. The corrections for this (allowing for the vertical profile of atmospheric 
density) require a knowledge of the spectrum of the emitted light. 

The uncertainty on the fluorescence yield remains the largest single contribu- 
tion to the overall uncertainty. In order to match the advancing capabilities 
of the UHECR detectors, it is being studied experimentally by several groups 
using different techniques |18|19ll20|21|22|23|24|25j . The yield and spectrum, 
as a function of atmospheric pressure, has been reported at several, often quite 
low, electron energies. A study of the light yield as a function of depth in elec- 
tromagnetic showers, and the sensitivity of the spectrum to depth, has been 
reported [26]. The present paper contributes light yield measurements over 
the range of pressures important for UHECR showers, and a survey of the 
spectrum. It makes use of a detector calibration technique systematically dif- 
ferent from other approaches. The goal was to reduce systematic uncertainties 
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in the fluorescence yield and spectrum below 10%, commensurate with other 
UHECR experimental uncertainties. 



2 Experimental method 

Much of the work reported so far has used radioactive sources to excite the 
air in a test cell. This corresponds to the low end of the dominant part of the 
shower's electron spectrum. The use of pulsed high energy electrons entails 
a different set of systematic issues, and has some advantages. For example, 
the monochromatic electron trajectories are easy to model, together with the 
flducial length for light emission from the test gas. With a pulsed beam, the 
light signals can be strong, statistics may be collected relatively quickly, and 
photomultiplier tube random dark noise does not contribute a background. 
On the other hand, linearity of the signal response must be checked, heavy 
shielding is necessary from stray radiation, and backgrounds must be studied. 

The beam available for this study was in the Final Focus Test Beam (FFTB) 
facility at the Stanford Linear Accelerator Center. Electrons at 28.5 GeV were 
delivered at 10 Hz in pulses 3 ps long. The apparatus was installed in an air 
gap in the beam vacuum line, with 50 micron thick stainless steel vacuum 
windows upstream and downstream. The beam particle trajectories were ef- 
fectively parallel, and their transverse distribution was measured nearby using 
transition radiation in visible wavelengths emitted from a titanium foil in the 
beam. This light was imaged by a CCD camera and image capture system. 
Beam spot widths were typically ~1 mm. 

The value deduced for the fluorescence efficiency depends directly on the mea- 
surement of the beam intensity. For this purpose, a toroid was mounted in the 
beam line upstream of the fluorescence apparatus [27]. The electron bunches 
passed nearly centrally through this ferrite ring which had an evenly spaced 
copper winding. The winding was coupled to front-end electronics, close to 
the beam line, which amplifled the current impulse and used a bandpass fllter 
to improve noise rejection. This signal was sent to the remote data collection 
system where it was digitized on every beam pulse. Beam intensities as low as 
10^ electrons per pulse could be measured with adequate resolution. At low in- 
tensity the pulse-to-pulse intensity variation could be ±30%. The calibration 
of this unit, including its linearity, was studied by injecting known charges 
into the electronics, and also, after removal from the beam line, by inductive 
coupling into the toroid itself of fast pulses simulating beam bunches. The 
calibration factor has been established with an overall uncertainty of better 
than 2.7%. Since the toroid was upstream of the fluorescence apparatus, it 
did not sense the perturbation on the electron flux caused by the thin beam 
windows. This has been evaluated by simulation codes, as discussed below. 
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The apparatus in which the fluorescence occurred, and was measured, is il- 
lustrated in Fig. [U It consisted of a 25 cm long, 15 cm diameter cylindrical 
stainless steel vessel, mounted coaxially with the electron beam. The upstream 
and downstream ends were closed off by 25 micron thick aluminum pressure 
windows. Inside this, a pair of 1.6 cm diameter thin aluminum cylinders were 
placed, one upstream, one downstream, coaxially with the beam, with a gap 
of 1.67 cm between them in the center of the apparatus. This gap formed the 
defined length of gas to be observed by the photomultiplier tubes. The internal 
surfaces were coated black to suppress scattered light, including the forward 
emitted Cherenkov radiation from the beam. 

Light from the gap could pass down two light channels which extended out 
radially through the cylinder walls. The channels were at right angles to each 
other, and were designated North and South. Their interiors also were black, 
and had baffles to suppress unwanted light paths. Each terminated, at 45 
cm from the beam axis, in a 1.2 cm diameter fused silica pressure window 
that formed the effective iris, as well as being the mechanical limit of the low 
pressure volume. The light continued on for 15 cm, undergoing a right angle 
reflection at a UV enhanced aluminum coated mirror. This turn in the optical 
axis allowed lead shielding to be placed between the photon detector and the 
beam line, eliminating the relatively direct radial path for scattered electrons 
or gamma rays to enter the face of the detector. 

A wheel of optical fllters was installed immediately after the right angle re- 
flection. Various narrow band fllters were available, as was a sample of the 
300 - 400 nm fllter used for the HiRes telescope, a clear aperture, and a blank 
position to study "dark" backgrounds. The rotation of the fllter wheels was 
controlled remotely. 

The flltered light was collected by one photomultiplier tube, (38 mm diameter 
Photonis XP3062 [28]) in each light channel. The PMTs were samples from 
the inventory of tubes used for the HiRes experiment. Their high voltages 
were monitored throughout the experiment. Diametrically opposite each of 
these light channels was a shorter cylinder in which was mounted an ultravio- 
let LED that was flashed between beam pulses to monitor PMT gain stability. 
Other LEDs were placed (outside the fluorescence light path) between the 
fllters and the PMT face. Along with the PMTs, in the same shielding enclo- 
sure, was placed a similar tube with a hood over its photocathode. This was 
intended to monitor noise not associated with fluorescent light, particularly 
the effects of penetrating radiation. Although the apparatus, and in particular 
the PMT section, was encased to the extent possible in lead shielding, some ra- 
diation could penetrate and excite the PMTs, depending on beam conditions. 
An additional hooded tube with a different gain was enclosed in a separate 
shielded cavity. 
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A gas system, with its controls outside the beam radiation enclosure, allowed 
the fluorescence cylinder to be fllled with dry air. For some measurements, 
ambient moist air was used to investigate the effect of water vapor. For sys- 
tematics studies, nitrogen or ethylene (which fluoresces only very weakly in 
the relevant wavelength range) could be introduced. The system pressure was 
varied in steps covering the range from 10 torr to 750 torr. 

The PMT and toroid signals were recorded using a standard CAMAC gated 
analog-to-digital-converter (ADC) system controlled by a PC. Also recorded 
were gas pressure, fllter position, PMT high voltages, and temperatures and 
gas humidity. During data collection, occasional triggers were imposed to mea- 
sure ADC pedestals and to pulse the LEDs used to monitor the PMT gains. 



3 Optical calibration of the detector system 

Various factors contribute to the light detection efficiency: the geometrical 
optical acceptance; the reflectivity and transmission of the optical materials; 
the PMT quantum efficiency; its amplification (gain); and the calibration of 
cabling and the charge-to-digital conversion efficiency of the ADCs. All the 
factors related to photons are wavelength dependent, although in the range of 
interest, the optical reflectors and windows are almost negligibly so. In order 
to measure all these factors in the intact optical set up, the system was moved 
to a laboratory where the vacuum windows were replaced by coated optical 
windows. Where the electron beam had been, the narrow beam from a nitrogen 
laser was aligned. The signal detected by the PMTs was, in this case, the 337 
nm laser photons that were Rayleigh-scattered from gas molecules in the 1.67 
cm target gap. Since only ~ 10~^ of the photons scattered in the target gap at 1 
atm., it was necessary to take great care to suppress stray light from elsewhere 
in the optical path. The laser beam was monitored by a probe calibrated by 
the manufacturer [29] to ±5%, and signals from the PMTs were measured 
with the system previously used at the electron beam line. By comparing the 
PMT responses to the LED flashers in the laboratory with those recorded in 
the electron beam enclosure, an average repeatability difference of 0.5% was 
measured, with a range from +2.5% to -1.4% for the different LEDs. For this 
contribution to the calibration uncertainty, we conservatively take ±2.5%. 

Measurements were made at a wide range of gas pressures between atmo- 
spheric and vacuum. It was found that the signal strength, normalized to the 
laser intensity, rose linearly with pressure. Fig. as expected from Rayleigh 
scattering. The intercept at the vacuum setting corresponded to the back- 
ground from errant laser rays, and the detected PMT signal from Rayleigh 
scattering was obtained from the slope. The actual Rayleigh scattering yield 
was calculated from the laser intensity and the gas density using the expres- 
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sions provided by Bucholtz [30]. A calculation based on Bodhaine et al. [31] 
gave results consistent within 0.2%. 

Additional uncertainties taken into account include those of environmental 
changes on the gas density (1.1%), and of the straight line fit to the pressure 
curve (0.2%). Measurements were also made with the HiRes wide-band filter in 
place. The ratio of signals with : without filter gave a value for its transmission 
at the 337 nm laser wavelength. This was also measured using a spectropho- 
tometer, and the discrepancy between the results of the two methods, ±1.8%, 
may be taken as an indication of the repeatability uncertainty. 

The calibration has been extended from 337 nm to the full wavelength range 
by comparing the DC current responses of the PMTs and NIST wavelength- 
calibrated photodiodes [32] , wavelength by wavelength, using the output from 
a monochromator. Fig. [31 Between 337 and 420 nm the PMT response was 
found to vary by less than 10%. Below 337 nm, however, it dropped off, as 
expected for a bi-alkali photocathode. By 300 nm it was as low as 30% of the 
337 nm performance. From wavelength to wavelength the rms fluctuations in 
the response measurement were always less than 2%. 

In measurements without optical filters, a few percent more light may be de- 
tected beyond the limits of the monochromator survey, and so the measured 
response function has been extended using manufacturer's data. The overall 
sensitivity to the air fluorescence spectrum has been computed using two ex- 
amples of measured spectra. One was that reported by the Airfiy collaboration 
[23] . and the second, from this experiment, is discussed below. In both cases, 
the spectrum was extended over the weak emission range to 600 nm, using 
the visible wavelength results from Davidson and O'Neil [T^. The results of 
the average response calculation for the two spectra agree within 0.1% where 
a HiRes filter is used, and within 3% for the case with no filter. Relative to 
the 337 nm light, for air fluorescence the system average response is 0.648 
when the HiRes filter is used, and 0.995 without a filter, and we assign an 
uncertainty of ±0.015. To explore the effect of a possible shift in the calibra- 
tion of the monitor photodiodes across the width of the spectrum, a factor 
proportional to the wavelength distance from the 337 nm calibration point 
has been introduced. By 420 nm the imposed change in sensitivity was 5%. 
Since most of the light is from wavelengths relatively close to 337 nm, however, 
this changed the calculated response to the fluorescence spectra by only 1.2%, 
which is taken as a contribution to the overall uncertainty. 

A correction was applied for the difference in geometry between the source 
of the light from Rayleigh scattering and that from the more wide-spread 
energy deposition from the electron beam. In the electron beam case, one 
source of loss is related to the limited transverse optical acceptance of the light 
channels. Another effect is the suppression of low energy scattered particles 
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at larger radii, caused by the material of the beam tubes and the rest of the 
detector, an effect that does not occur in atmospheric showers. The spread 
of the energy deposited by the electron beam was simulated in a model of 
the apparatus using EGS4 [33j, and the acceptance efficiencies at the optical 
iris were calculated numerically for both the laser and electron beam cases. 
The factor obtained from Rayleigh scattering for the efficiency of converting 
photons to ADC counts must be reduced by (3.2±0.25)% for the case of the 
electron beam in the apparatus fiducial volume. In addition, the ratio between 
the energy deposited in a 1 cm length of free air without beam windows and 
that in the fiducial volume is 1.0837 ib 0.0015. For calculations using the 
actual energy deposit in the fiducial volume, an overall simulation systematic 
uncertainty of 1% should also be applied. 



4 Data analysis 

4-1 Data processing and background subtraction 

The data were accumulated in runs of several thousand beam pulses, and 
the gas pressure or optical filters were changed between runs. Within each 
run, the ADC signals were corrected for the zero-beam pedestals. The signals 
from the active PMT and the background counters were plotted, pulse by 
pulse, against the measured beam charge, as in Fig. HI The slope of the PMT 
versus beam regression is proportional to the fluorescent excitation caused by 
the electron beam, together with a beam-related background contributed by 
radiation penetrating the shielding and exciting the PMT. 

For this work, only the results of the north PMT will be discussed. The south 
counter, with a higher gain selected for use with narrow wavelength-band opti- 
cal filters, was in saturation for the configurations with relatively transparent 
filters used for overall sensitivity determination. The narrow band filter results 
are intended to be the subject of a future report. 

At the low beam intensities used, the delivered pulse charges typically were 
variable by ~ ±30%. This permitted the signal vs. beam plot to be studied for 
systematic effects. For a selection of runs with a suitable low beam intensity, 
where the PMTs were checked to be responding linearly, the regression was 
found to intersect the beam axis at (—0.41 ± 0.31) x 10'' electrons per pulse, 
consistent with zero. As a consequence, the origin has been applied as a con- 
straint on all subsequent fits, where the typical beam intensity was ~ 10^. To 
allow for the uncertainty, a conservative contribution of 1% has been included 
in the overall uncertainty evaluation. 
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Examination of the data allowed the limits of linearity to be determined. Two 
effects have been observed. The normal "saturation" sensitivity fall off of the 
PMT was seen to affect a small fraction of high pressure data where no op- 
tical filter was used. This was avoided by imposing a cut on beam intensity 
in runs showing this effect. A second, and more important, effect has its ori- 
gin in non-linear effects caused by the collective electric field impulse of the 
intense electron bunch. This field has the ability to accelerate electrons from 
ionization events, and these in turn cause molecular excitations, adding to the 
fluorescence signal. Related effects have been studied in the same beam line 
|36] . This enhancement of the fluorescence is illustrated in Fig. O We find 
that the beam pulse intensity threshold for the enhancement decreases with 
gas pressure, from 1.5 x 10^ at 1 atm to 0.8 x 10^ at 50 torr, for beam spots 
X (Ty of 1 mm X 1 mm, and 1 mm long. Beam intensity cuts were made to 
avoid this problem, and, as a result, some higher intensity runs could not be 
used. 

The backgrounds caused by radiation penetrating the shielding were mea- 
sured directly. This was achieved by rotating the filter wheel to the opaque 
position and fitting against beam intensity as before. Runs in this condition 
were interspersed among the others. Short term variations in the background 
were accounted for by monitoring the signals of the two permanently hooded 
"background" PMTs. Although their pulse-to-pulse fluctuations were large, 
over run-length intervals their averages tracked each other well. They have 
been averaged to provide a run-to-run adjustment of the opaque-filter back- 
ground runs. These backgrounds, typically ~ 5% for data without a filter, 
or ~ 7% with the HiRes filter, were then subtracted from the results of the 
neighboring fluorescence run fits. 

There were electron beam runs in both 2003 and 2004, and the Rayleigh scat- 
tering calibration was performed after the latter. The stability of the system 
between the two beam runs was tested by comparing sets of data taken under 
similar conditions, and the variations observed lead us to assign an uncertainty 
of ±2% for the long term stability. 



-^.^ Photon yield in dry air 

With any given setting of pressure and filter, several runs were taken, sepa- 
rated in time by hours to days. At each setting, the results show very good 
consistency, and so they have been averaged. The variation between them has 
been used to estimate a run-to-run uncertainty, which is included in the final 
results. 

After applying the calibration discussed above in the optical calibration sec- 
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tion, the dry air results are given in Tabled] for some pressures in the range of 
interest for UHECR detectors. Only HiRes filter data are available at 50 torr. 

The uncertainties that should be applied to these figures are listed in Table [2] 
in units of percent. They have been discussed in the narrative above. 

From these, the overall uncertainty on the yield per MeV is 7.5%: i.e. at 760 
torr the yield is 20.8 ± 1.6 photons per MeV. 

Results have frequently been quoted in terms of photons per meter, the elec- 
tron energy loss per meter being a function of its energy. We also give these 
values in Table [3] for both no filter and HiRes filter cases. The HiRes filter 
results are also shown in Fig. [61 For air, the fluorescent yield of an electron 
track barely increases with increasing pressure above ~ 0.1 atm. This is be- 
cause of increasing competition from molecular collision processes, principally 
involving oxygen, which compensate for the increased energy deposit with gas 
density. The fit line shown is the expression aP/{l + bP) P^1[2U] motivated 
by quenching of the fluorescent molecules by pressure dependent collisions. 
Here P is pressure and a, b are constants. The error bars shown in this case 
are the point-to-point relative uncertainties of 1.42%. The overall 7.5% scale 
uncertainty still applies. 

The data were taken at a temperature of 304 K. The effect of tempera- 
ture T is to change the collision rate with the quenching oxygen molecules, 
represented as a modification in the denominator of the above expression: 
aP/(l + 6P(T/304)i/2) [X3]l20] . The increase in the signal would, for example, 
be 1.3% at the 296 K temperature reported in Ref [22]. At 273 K it would be 
5.4% , and would reach 10% only at 250 K. 

4-3 Nitrogen, ethylene and humid air 

As a systematic check, some data were collected with nitrogen in place of air. 
CoUisional de-excitation is much weaker in pure nitrogen than in air. In the 
absence of heavy optical filtering, the stronger nitrogen signal at higher pres- 
sures caused non-linearity in the PMT at all but the weakest beam intensities. 
In response, we have fitted a quadratic function to the nitrogen data. The fits 
were repeated for various beam intensity ranges up to limits between 1.2 and 
1.5 xlO^, and the variation in the linear term was noted. This variation was 
taken as an indication of the systematic error of this approach, 4% for data 
without an optical filter, 0.2% with the HiRes filter. Additional uncertainties 
came from evidence of a small manifold leak of dry air into the nitrogen (1.7%) 
and from the signal attenuator (4%). It has been found that, at 750 torr, the 
signal ratio between nitrogen and dry air is 6.51 ± 0.37 with no optical filters, 
or 6.84 ± 0.29 with the HiRes filter. This is in agreement with our previously 
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reported [21] value of 6.6 ± 0.2. 



Results from a short run with a fill of ethylene were processed in the standard 
way, but the signals were barely above background. At 750 torr, the strengths 
relative to air were (0.34 ±0.13)% with the HiRes filter, and somewhat larger, 
(0.88 ± 0.09)% without a filter. This can be used to set an observational limit 
on the Cherenkov light background. 

Ethylene yields 2.5 times as much Cherenkov light as air. If all the ethylene 
signal were from this process, then its contribution to the air signal, at 750 
torr and with the HiRes filter, would be (0.138 ± 0.052)%, or < 0.21% at 90% 
confidence. The wider wavelength range of the no-filter case would not increase 
the Cherenkov signal by more than xl.4 relative to the HiRes filter case. This 
is already negligible for our purposes, but it is likely that much of the weak 
ethylene signal is the result of fluorescence. From its design parameters, the 
black beam tube is expected to suppress the Cherenkov light contribution at 
1 atm, the worst case, to < 0.1% of the air signal. 

Filtered air from outside the building, containing a relatively large fraction 
of water molecules, was also studied. At the overall pressure of 750 torr, the 
partial vapor pressure of water was 11 torr, or 1.5%. The data taking and 
processing were similar to those for dry air data. The results show that, at 
750 torr the humidity reduced the light yield by (7.4 ± 0.3)%, consistent be- 
tween HiRes- and no-filter cases. Data from a single no-filter run and a HiRes 
filter run, with room air at 245 torr, were also available. In these cases the 
yield reductions were (3.5 ± 1.4)% and (8.8 ± 2.6)%. We do not consider the 
discrepancies to be significant, and so take a weighted average of both pres- 
sures. That is, 1.5% water vapor content (equivalent to 7.1% of the oxygen 
partial pressure) suppressed the signal by (7.3 ±0.3)%. By comparison, at the 
temperature of the standard atmosphere at 400 torr or 5000 meters altitude, 
the water vapor partial pressure saturates at 0.25%. 



5 Spectrographic observations 

During part of the data taking, an opportunity arose to use an independent 
arrangement in the beam line to record emission spectra with relatively high 
resolution. This was accomplished by mounting an 8 cm long, thin-windowed, 
gas cell in the beam, and observing the fluorescence light emitted perpendicu- 
larly to the beam through an aperture 2.6 cm along and 1.3 cm transverse to 
the beam. The gas cell was baffled to suppress scattered Cherenkov light. Two 
aluminum-coated 45° plane mirrors, a fused silica pressure window, and a fo- 
cusing mirror were used to match the light from the cell to the f/4 acceptance 
of the 120 mm focal length spectrograph [51] in a heavily shielded enclosure. 
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Signals were measured pulse by pulse by using a multi-anode photomultiplier 
tube with a linear array of 32 pixels spaced at 1 mm [35]. The anode signals 
were digitized by the standard CAMAC 11-bit ADCs. Beam intensity infor- 
mation was not available synchronously for this data set, and so the data does 
not support comparisons of total yield at different times, although the shapes 
of spectra may be compared. 

Between runs the light path was deflected away from the spectrograph slit, and 
this gave a background measurement for off-line subtraction. The raw signals 
were averaged over each run taken at a fixed pressure, and the background, 
including the ADC pedestal, was subtracted. The plot of these values already, 
at this raw stage, showed the characteristic spectrum of the gases studied, 
dry air (including argon and carbon dioxide) and nitrogen. Unpopulated gaps 
between the expected spectral lines excluded the presence of a significant 
continuum spectrum such as a Cherenkov radiation background. 

After the beam run, the wavelength calibration of the system was made by 
measuring the positions of the lines from a mercury discharge lamp, and the 
sensitivity as a function of wavelength was studied by recording its perfor- 
mance when illuminated by the standard continuous spectrum of a deuterium 
lamp. It was found that the response of the system as a function of wavelength 
was dependent on the light path through the spectrograph. When using a cal- 
ibrated light source, it was not possible to match the distribution in space 
and angle of the beam- induced fluorescence light, and this has led to an un- 
certainty of ~ 15% in the relative light yield of the short wavelengths ~ 315 
nm, relative to the rest of the spectrum which extended to 415 nm. 

The analysis corrected for crosstalk between neighboring anodes, given as 3% 
by the manufacturer, consistent with our observations. Additionally, there 
was evidence for a loss of linearity of anode signals for the strong lines, also 
consistent with the manufacturer's specifications. A correction for this was 
evaluated using pulse amplitude data for the strong lines at 337 nm and 357 
nm and the line at 391 nm, using different pressures and a wide range of 
approximately known average beam intensities. Most of the significant lines in 
our spectral range are from neutral nitrogen excitations of the second positive 
(2P) system C^Hu — ^ B^^g- The non-linearity correction for the 2P lines was 
found to rise with pressure to a plateau above ~ 300 torr. This was expected, 
since the light flash width has been observed to become narrower inversely with 
pressure |2T], thereby increasing the peak dynode current. The 391.4 nm band, 
labeled 1N(0,0), is a transition between vibrational levels i?^E+ and X^S"*" of 
ionized nitrogen. For the IN line, the pulse width had been observed to have 
already become narrow by 50 torr, and, as this led us to expect, the non- 
linearity correction rose more quickly with pressure to a constant level. For 
stronger lines at higher pressures, the non-linearity corrections were ~ 10%. 
An average systematic uncertainty of 17% of the correction was assigned. 
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An example of a spectrum, from air at 155 torr, corrected in this way is given 
in Fig. [71 Each bin corresponds to one anode and the sum of the bins is unity. 
In Fig. [S] are shown, for air at various pressures, the relative signal strengths 
within several wavelength ranges incorporating the most important lines. At 
each pressure, the corrected signals are expressed as a fraction of the total 
signal at that pressure. There is not a strong change in the spectrum over 
the pressure range of interest (above ~100 torr). The transmission efficiencies 
of these spectra through a HiRes filter vary by 0.46% rms, the largest devia- 
tion from the mean being 0.70%. Also, the Rayleigh scattering transmission 
efficiencies through 30 km of 1 atm. clear air varies by less than 1% for the 
spectra at pressures above 60 torr. Associated with its increased decay time, 
the 391 nm transition of the IN system is seen to become relatively more 
important as pressure decreases below 60 torr. This is expected theoretically 
[T5] because of differences in the radiative and coUisional de-excitation rates 
of this quantum state. 

In order to translate from a "spectrum" of anode signals to one of emission 
lines, an additional step was required. This was because the sensitivity profile 
across the 1 mm width of each anode rose and fell between the inter-anode 
gaps. It followed an elliptical shape, averaging a 20% loss from uniform sensi- 
tivity, while dispersed images of the slit could fall into poor acceptance regions, 
increasing uncertainties. A parameterization of the manufacturer's average ac- 
ceptance profile for this was included in a numerical model of the device. This 
permitted the alignment and slit width calibrations to be applied and varied 
within the estimated uncertainties. 

The air fiuorescence spectra for resolved emission lines between 300 and 415 
nm were extracted from this simulation. The result at 155 torr, approaching 
the pressure of the upper altitude range of interest to UHECR detection, is 
illustrated in Fig. [HI where the signals are expressed as a fraction of the total 
yield in the full wavelength range. At this pressure, the photomultiplier tube 
non-linearity corrections, and their uncertainties, are relatively small. 

It may be of significance that the optical system focused a region of just ±2 
mm around the electron beam on to the spectrograph slit, and so excluded 
fluorescence from energy deposited beyond this. It appears possible that the 
spectra from the core region of the beam and the peripheries may be somewhat 
different p!6]. The 391.4 nm band, labeled 1N(0,0), is caused dominantly by 
direct ionization by the high energy beam, while the 2P system, responsible for 
the other important lines, is excited by low energy scattered electrons. With 
the tightly selected field of view, the 391 nm region may have been detected 
more efficiently than the rest of the range. 

In data runs taken with a lower resolution grating, a comparison of dry air and 
room air, with 1.5% partial pressure of water, showed no discernable humidity 
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related change in the spectrum profile. 



6 Discussion 

Several points are worthy of attention. The sensitivity calibration of the pho- 
tomultiplier system with the HiRes filter has an uncertainty better by a factor 
of ~ 2 than has been available before. The change in the HiRes UHECR 
spectrum caused by this is within previous uncertainties, as will be reported 
elsewhere. 

It is also reassuring that the observed spectral structure of the emitted light in- 
deed comes from the expected nitrogen bands [T3p!7] . In comparing with spec- 
trum measurements reported elsewhere, one finds general agreement. There 
are discrepancies in detail, however, reflecting the difficulty and range of pa- 
rameters and techniques involved. In Fig. [10] are illustrated various reported 
signals in wavelength bands suitable for the UHECR detector range 300 to 
420 nm. Reports from the early work of Bunner [13], the optical filter work 
of Nagano and collaborators [12], the Airfly collaboration [23] and from this 
paper, are shown for comparison in the same bands. The values plotted have 
been expressed as a fraction of the total signal reported within the range. 

The agreement is adequate for present UHECR data analysis. The transmis- 
sion efficiency of the HiRes optical filter for these spectra is the same within 
1%. We have compared the effect on them of the wavelength dependence of 
Rayleigh scattering. The transmission values for these four spectra, calculated 
through a 10 km air column at 1 atm., average 0.481, with a worst case de- 
viation from the mean of 0.009. At 20 (30) km atm the average transmission 
is 0.242 (0.124) with a worst case deviation of 0.008 (0.009). (The average 
pressure along a light ray from the peak of a shower to the UHECR detector 
is typically 0.6 - 0.7 atm.) At 30 km atm, these differences, amounting to up to 
7% of the transmitted signal, will become of significance for future detectors, 
and will then need clarification. 
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Table 1 

Photons per MeV as measured with no filter and penetrating the HiRes-filter. An 
overall uncertainty of 7.5% applies to the yields (see text) 

photons / MeV 



pressure (torr) no filter through HiRes filter 



760 


20.8 


14.0 


495 


32.0 


21.7 


242 


64.3 


43.0 


97 


157.6 


105.2 


50 




182.2 
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Table 2 

Contributions to uncertainty on the photon yield, in units of percent 



uncertainty contribution 


% 


beam calib. 


2.7 


signal splitter 


1 


zero constraint of fits 


1 


run-to-run stability 


1 


laser vs e-beam light source shape 


0.4 


simulation 


1 


spectrum sensitivity, open filter 


1.5 


spectrum sensitivity, HiRes filter 


1 


beam line vs lab stability 


2.5 


2003 data calib. 


2 


filter consistency checks 


1.8 


PMT relative spectral response 


1.2 


Raylcigh scattering: 




laser power 


5 


gas density for laser scattering 


1.1 


theoretical calculations 


0.2 


fit slope 


0.2 
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Table 3 

Photons per meter for no filter and HiRes filter cases. The relative uncertainty 
between points at different pressures is 1.42%, and an overall scale uncertainty of 
7.5% applies to all points. 



pressure (torr) 


no filter 


HiRes filter 


750 


5.059 


3.413 


495 


5.029 


3.403 


242 


4.848 


3.240 


97 


4.686 


3.128 


50 




2.784 
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J3J1 




Fig. 1. Sectional view of the apparatus. The electron beam axis is indicated by the 
arrow. The optical path is perpendicular to the electron beam along a baffled tube, 
through a fused silica pressure window, to a right-angle reflection. Following this is 
a wheel of optical filters, LEDs (outside the fluorescence light path) for monitoring 
stability, and then the PMT. The opposite arm contains an on-axis LED. 
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4 4 

x10 x10 

~ pressure 

Fig. 2. PMT response against Rayleigh scattering intensity, controlled by changing 
air pressure. The ordinates are ADC counts, corrected for pedestal and normalized 
by the laser power measurement. 
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Fig. 3. PMT response relative to that at 337 nm, and transmission of HiRes filter, 
vs. wavelength. 
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Fig. 4. Fluorescence signal plotted pulse-by-pulse against beam intensity at 750 torr, 
HiRes filter 
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Fig. 5. Fluorescence plotted against beam intensity at 250 torr, no optical filter, 
showing the enhancement occurring at higher beam intensities. 
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Fig. 6. Photons per meter measured using tlic HiRes filter. The error bars shown 
indicate the point-to-point uncertainty of 1.42%. An overall scale uncertainty of 
7.5% applies to all points. See the text for the fit expression. 
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Fig. 7. Fluorescence spectrum of 155 torr air as seen by the multi-anode PMT. Cor- 
rections for transmission, conversion efficiency and non-linearity have been applied. 
The sum of entries is normalized to unity. 
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Fig. 8. Pressure variation of the relative contributions of various spectral ranges. 
The most noticeable structure is the rise below 60 torr of the fraction of light in the 
391 nm IN band, illustrated by the line. 
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Fig. 9. Line spectrum from 155 torr air, obtained by matching known lines with the 
response of the spectrograph. The sum of the hne strengths is set to unity. 
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Fig. 10. Comparison between published emission strengths and those reported here, 
grouped into six wavelength bands. Minor adjustments have been made to allow for 
compatible wavelength ranges. In each case the sum of emissions is normalized to 
unity. 
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